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I. INTRODUCTION 
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^ ■ The propagation of electromagnetic fields in curved space-time has always been an important problem since the 
formulation of general relativity. Gravitational Doppler shift, deflection of light rays, and gravitational lensing have 
been useful tests for general relativity and have represented important tools in astronomical observations (see for 
instance [0J|]). 

In the first approach to this problem, the gravitational background was considered as a sort of anisotropic medium 
^ ■ (a good introduction to this kind of approach is found in ||^ , while a recent exhaustive treatise is |^ ; for the passage 
■rj" I from wave optics to geometric optics in curved space-time see |^). In the literature, on the level of geometrical 
. optics, many authors pointed their efforts at gravitational fields of isolated physical systems (such as a star at rest or 
■ a rotating body) which were not radiative ones (see |^|ll[]). It was also considered the problem of light fluctuations in 
a gravitational wave background (see [ p2| ) ; however in this work the author did not address his attention to possible 
polarization effects. 

In the following years, in the framework of the linearized theory of general relativity, much work has been done 
^ to describe the propagation of electromagnetic fields in a gravitational wave background, using a different approach 
with respect to the usual one. In fact many authors have split the electromagnetic tensor (or the 4~vector potential) 
I in a sum of two terms: the first one is the flat space-time solution, while the second term describes the perturbation 
^jpj' due to the weak gravitational field (see ||l^-|20|). Nevertheless, as it will be shown in Sec. II, in this way the original 
1^ • . physical problem is changed. It is for this reason that we think that the former approach is the most suitable for 
the problem we are concerned with. An attempt in this direction has been recently done (see [^) but the Lorentz 
condition is not properly handled, being fulfilled only for particular directions of propagation and amplitudes of the 
5J] • 4-vector potential. 

Ci It is therefore seen that, as far as the problem of the propagation of electromagnetic radiation in the field of a 

gravitational wave is concerned, a satisfactory solution is lacking which takes into account the tensor nature of Maxwell 
equations, even in the geometrical optics limit. Purpose of this paper is to investigate this problem. Therefore, besides 
the usual phase perturbation, the consequences on the polarization of the electromagnetic field are evaluated. It is a 
result of the present paper that, when Lorentz condition is met, nontrivial results in the polarization perturbations 



are found (see Sec. VI and the conclusions). 

As it is well known. Maxwell equations are a set of redundant equations. If the problem is considered from the 
electromagnetic 4-potential point of view this redundance could be partially eliminated. The free electromagnetic 
equations in curved space time for the vector potential (de Rahm equations) read 

A'^-'^^ - A"./ = (1.1) 

(all through the paper, notations and conventions as in ^J). In order to get a simple set of equations, adopting the 



standard approach, we impose the Lorentz gauge condition on vector potential, bringing Eq. (1.1) into the system 

Ai"'",, = 



A^.„ = 



(1.2) 



Once a solution of the previous system is found, the antisymmetric tensor of the electromagnetic field, which is the 
physical measurable quantity, is obtained by the usual relation 



1 



A,. 



(1.3) 



From the above equations it is clear that the direct interaction of a gravitational wave with an electromagnetic field 
is due only to the metric tensor. This means that in curved space-time the differential equations of electromagnetism 
have variable coefficients. However the explicit form of the dependence from the metric tensor depends critically 
on the reference frame used to perform calculation. Our strategy is to consider system ( |1.2D in a reference frame 
in which calculation is easier. This system is the so called transverse-traceless (TT) gauge (see 0]). However the 
reference frame in which measures are performed is the Fermi Normal Coordinate (FNC) system (see for instance [Q). 
Therefore once the antisymmetric tensor of the electromagnetic field is known in the TT gauge we find its expression 
in FNC by means of the usual transformation rules (the electromagnetic field tensor is indicated by F^^, in FNC and 
by J-ap in TT gauge; are the FNC while y" are the TT gauge ones): 



F, 



dy" dy^ 



(1.4) 



The connection between the two frames is known for every point of the space-time in the linear approximation of 

weak gravitational fields. This result is found in p2| , p3| . 

The paper is organized as follows. In Sec. ^ we find a solution to the de Rahm equations (1.2) in TT gauge. In 
and Sec. IV we write the electromagnetic tensor components in TT gauge and FNC respectively. Finally in 
we apply our results to the calculation of the polarization effects. 



Sec. m 
Sec. VI 



II. SOLUTION OF THE DE RAHM EQUATIONS IN TT GAUGE 

For weak gravitational field, to which we shall limit ourselves in this paper, the metric tensor could be written as 

5m"^ = '7m'' + ^^M"^' << 1 (2.1) 

where rj^^ is the Minkowski metric tensor with positive signature. In this approximation the particular coordinate 
frame called TT gauge (see [Q) is characterized hy h — rj^'^h^j^y — 0, ft-o/j — ^ind h^-' j ~ 0. Neglecting second order 
terms in h^j^^ the linearized homogeneous de Rahm equations with the Lorentz condition are 

A'^-l - h'-^A'-^^^p + {h^J' + h^^,^ - h£^^^) = 

(2.2) 

This system is a set of differential equations with variable coefficients. In the approach of fl^-p^ the vector potential 
A^ is split into two parts: 

= + ^A^ (2.3) 

where the first term is a solution of the flat space-time wave equation while the second one represents the perturbation 
due to the weak gravitational fie ld. In this approach all the quadratic terms in the perturbation vector are neglected. 



With this assumptions system (2.2) becomes a set of differential equations for the perturbed vector potential with 
constant coefficients (the d'Alembert operator) and a source term which involves the unperturbed potential and /i^i/ 
(e.g. [p^l ). It is our opinion that this procedure could in general lead to some misunderstanding on the nature of 
the interaction when the "unperturbed" electromagnetic field is oscillatory (this is not the case of |l^ in which the 
"unperturbed" electromagnetic field is constant). To show this point let us consider the case of an ordinary differential 
equation, the Mathicu equation (sec for instance [^,^), which may be written as 

q + uj'^[l + h{t)]q^O (2.4) 

and which describes, for instance, the charge on the plates of a condenser with a variable capacitance in an oscillating 



LC circuit iQ. Supposing h{t) << 1, a splitting q = ^q+\ [in which is the "unperturbed" solution of Eq. (2.4) 
when h{t) — and \ is the perturbation due to the presence of h{t)] brings to the equation of a harmonic oscillator 
for \ with a driven force proportional to the product of h(t) and %. It is a matter of fact that, in general, the new 



set of equations is not equivalent to Eq. (2.4) ||2j|. The physical reason lies in the fact that splitting of the charge 
changes the original problem of a free circuit with variable capacitance in a different problem of a circuit with constant 
capacitance and an external electromotive force. In the same manner, splitting of the vector potential changes the 
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problem of the propagation of a free electromagnetic field in a gravitational background, acting as a sort of anisotropic 
medium, in the problem of the generation of a field by a given current. For these reasons our approach will be that 
of performing calculations without splitting of the vector potential. This strategy was yet carried out for the problem 
of the interaction of a gravitational wave with electromagnetic circuits in [^6| . 

Let us now state precisely the physical problem which we are intended to solve. Consider a TT gauge reference 
frame in which a plane gravitational wave propagates in the z direction. In this system the only non-vanishing 
components of /i^^ are [<i>g — x{y^ — y^)]- 

hn{<fg) = -h22{^g) = h+{<^g); /ll2(*g) = /l2l(*9) = /»x (2.5) 

Let us consider now a free electromagnetic field which in the absence of gravitational waves propagates in a given 
direction described by the spatial component of the four wave vector fc^. The vector potential for such a field is given 
by — C^^g{k^y^), where are four constants that satisfy Lorentz condition C^fc^ = 0. A gravitational wave, 
modifying the geometry of the space-time, acts as an anisotropic medium whose dielectric and magnetic properties 
change in space and in time, 1/x and l/(cx) being the order of magnitude of space distances and time intervals for 
which variations take place appreciably. If 

|fco|=fc»X, (2.6) 

then the amplitude, polarization and direction of the electromagnetic field remain practically constant over distances 
of the order of 1/fc and times of the order of l/{ck). Therefore we can assume that the form of a solution of Eqs. ( |2.2| ) 
would be 



(2.7) 



where a'^ is a 4-vector function of space and time and X is a scalar function. In the framework of assumption ( |2.6| ) 
we can suppose that first derivatives of and second derivatives of I are small quantities with respect to and 
fc^ respectively, so we can neglect superior order derivatives. This approach is one order of approximation beyond 
the geometrical optics limit allowing therefore to find not only the phase shift (that is to say the scalar function X, 
the eikonalL but also the change in the polarization of the electromagnetic wave (for all these considerations see for 
instance [p7[). 

In order to solve the problem it is important to find, within the Lorentz condition, an electromagnetic gauge for 
which the equations are eas ier t o solve. One can see immediat ely that, because of the space-time dependence of the 
gravitational wave [see Eq. (2^)], the first equation of system (2^) is the same for ^ = 3,0. Therefore we fix a gauge 
for which 



^3 = A" 



(2.8) 



With this choice the equations for A^' (from now on indices r, s, t run from 1 to 2) are independent on — A'^ . Once 
A^ are found, by means of the Lorentz condition one can find immediately A'^ = A° which fulfill also the first equation 
of (2^). To carrying on this program it is convenient to perform a coordinate change. We put 



T = 


iy' 


-2/°)/V2 


X = 


y' 




Y = 


y' 




Z = 


iy' 


+ 2/°)/V2 



(2.9) 



The new metric tensor G^i/ can be written as [see Eq. (|l.4[)] 



G 



QyO- Qyf^ 



(2.10) 



where 



N.,, 



/O 1 

10 

10 

V 1 , 



(2.11) 



If we call A'^ and K^^ the vector potential and the four wave vector in this new frame we get 
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A'' = 



( (A3 - A")/V2 \ 
\ (A3 + A")/^/2 / 



V2 



V2 



(2.12) 



In this frame the fixed gauge is written as the usual Coulomb gauge in flat space-time, that is to say: 

dA' 



A"^Q; 



dX' 



= 0. 



The first equation of (2.2) for the t components is written as 



a, 13 



while the equation for A^ is given by 



h\flA\^ = 



Once A'' are found, from Eqs. ( 2.13 ) we can find a general expression for A ; we can write 



A-" 



dA^ 



dZ' + g{X,Y,T) 



where, in order to fulfill Eq. ( 2.15D , g{X,Y,T) must be a solution of 



d^g{X,Y,T) „ d^A'' 



dX^dX' 



= 2 



dX^ dT 



+ V.o(V2xT) 



dA"^ 



dX^ 



(2.13) 
(2.14) 
(2.15) 

(2.16) 
(2.17) 



z=o 



As we shall see later this function could be chosen in such a manner to have a solution which describes "plane waves" 

travelling in some direction. 

In order to find a solution of Eq. (^.14 ) we exploit consequences of assumptions ([2 .61) and (2.7). For what was said 
we neglect second order derivatives of a'', and terms like h^^,a<i^ ^. Therefore Eq. (2.14) is divided into the following 
equations for and X: 











(2.18) 
(2.19) 



where /ir? means the n-order derivative of h^s with respect to the argument \/2xT . 

Now we notice that the "unperturbed" polarization amplitude and eikonal of the electromagnetic field are not 
oscillating quantity: indeed the polarization amplitude is a constant and the eikonal is a big quantity, changing by 
27r when we move through one wavelength, while the gravitational perturbation occurs on the scale of \/x » ^/k. 
Because of the smallness of the metric perturbation [see Eq. (2.1)] and of its coordinate dependence we can therefore 
safely put 



= S + u{T) 
+ b\T) 



(2.20) 
(2.21) 



where S — k^x"' + is the usual phase of the flat space-time, are constant and finally u{T) and 6*(T) are the 
small perturbations due to the gravitational wave background. With these assumptions Eq. (2.19) become 



1 KrK, 

2 " 

where dot means derivation with respect to T . The solutions to these equations are simply given by 

1 KrKs 



i/'-^(V2xT)-i/'-^($o) 



uo 



(2.22) 
(2.23) 

(2.24) 
(2.25) 
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where H'^^{x) is a primitive of h'^^{x). As far as is concerned, from Eqs. (2.16), (2.20), and (2.21) we can write 



J<3 



(2.26) 



where S'o = Sl^^^. Let us now fucus our at tentio n on the function g{X, Y, T). Taking into account Eqs. ( ^.22 ), (2.23), 
and the relation N""^ K^Kp = 0, then Eq. (2.17) could be written as 



52 



92 



A solution of this equation is 



,_d d_ 



^3 



g[T) = -^^[B^ + ¥{T)]f{S,+u) 



Putting this solution in Eq. (2.26) we find 



Kr 



A' = --f[B^ + b^{T)]f[S + u{T)] 



(2.27) 



(2.28) 



(2.29) 



The solution described by Eqs. ( 2.13| ), ( ^.20 )-( 2.25| ), and ( ^.29 ) represents a "plane wave like" solution for the 
electromagnetic field. 

Getting back to the TT gauge system in which we wrote Eqs. (v^) we finally find 



f{S + u) 



(2.30) 



A^ = A' = - 



ks + ko 



B^--h^,i%)B^+b^o 



fiS + u) 



where 



1 k g 

It = — - 



2x(fc3 + fco) 



(2.31) 



is the phase shift function. We have therefore written a plan e wave like solutions of Eqs. (2.2) in the particular Lorentz 
gauge for which A^^ = A'^. As one can see from Eqs. (2.30) the amplitude of A^ = A" tends to infinity when ko + 
goes to zero (that is to say for parallel interaction). However this is only a result due to the particular chosen gauge, 
without any physical meaning. In fact in the next section we shall see that the electromagnetic field does not present 
such infinity problems. 



III. ELECTROMAGNETIC FIELD IN TT GAUGE 



In this section we find the components of the antisymmetric tensor of t he electromagnetic field in interaction with 
a gravitational wave in the TT gauge. To this aim we shall use Eq. (1.3). In order to do that we must before write 
the covariant component of A^^. We have 



A3 = A3 = AO = -Ao 



(3.1) 



We therefore have: 



A3 = -Ao = -T^iB^ ~ \h\B^ + &S) f{S + u) 



(3.2) 



Ar = {B^ + -h^,B^ + hl) f{S + u). 
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In order to be consistent with the approximations made [see Eq. (^^) and its consequences] we neglect first order 
derivatives of /i^^ in Eq. (O). The components of the electromagnetic field tensor in the TT gauge are therefore 
written as 



and 



Hi = -kr ( B'' - ^h\B' + % ] f'{S + u) 



+ 



fco + fcs 



B^~-h%B'+¥, 



fiS + u) 



k, B"- + -h\B' + 6^; + u.B 



^12 
•?>3 



where we have used the fact that u^q = — and 

",3 = 



kAB^ + -h\B^ + hl\-k2{B' + -h\B^ + hi 



fiS + u) 



fco + fcs 



B^--h\B*+b^, 



fiS + u) 



2fco + fc3 ^ 



(3.3) 



(3.4) 



(3.5) 



IV. ELECTROMAGNETIC FIELD IN FNC SYSTEM 



The aim of this section is to write the expression for the electromagnetic field tensor in FNC. In the framework of 
linear approximation, TT-coordinates are given by FNC by means of 



(4.1) 



which involves the infinitesimal vector field e^; this vector field is known in every point of the space (see |2^) therefore 
solving the problem of the connection between FNC and TT-gauge. If the TT-gauge perturbation of the Minkowskian 
metric is a plane wave propagating in the z direction, with angular frequency x then the infinitesimal vector field e'^ 
may be written as (see p3[|) 



1 

2^ 



Xx^ 



XX'' 



- ^hrshxx")^ X" 



(4.2) 



where v — 0,3, and we remind that Hrs{4') is a primitive of hrs{(j)) and that r, s = 1, 2. Obviously a plane wave is only 
an approximation of the gravitational field in a relatively small region around a point far away from the source. This 
implies that (^|^) represents the infinitesimal vector field only in the region in which the plane wave propagation 
is valid. In which follows we consider source distances and gravitational wavelengths for which Eqs. (4.2) hold in the 
region of interest. 

The electromagnetic tensor field components in FNC are given in term s of the TT gauge ones by means of Eq. (1.4). 
If the transformation rule between the two system is given by Eq. (4.1) then we have 



(4.3) 



In the range of validity of Eqs. (4.2) for plane monochromatic waves from the z direction we can write [a = x(x 



-XX" 
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Fat — J-Q 



7z 2^' 



Hst{a) - Hst{ao) 



Xz 



- hstiao) 



03 



2x3 



hrs{a) — hrs{ao) ^f^^^ 



XZ 



(ao) 



V(a)-V(ao)_l^(,) ' 
XZ 2 



f 

^ St 



(4.4) 



^03 — ^03 H 



ff/(a)-F/(ao) 



V(a) 



Hrs{a) — Hrs{ao) 



J'Os 

,3 



- Ks{ao) ^h^^isM 



03 



X^ 



s3 



(4.5) 



F\2 — P\i H — 



s3 



g/(a)-g/(Qo) 



V(ao) 



(4.6) 



Ffi — 

x"" 



-z 2^'* 



r3 



ij'rt(a) - Hrtjao) 
XZ 

H/{a)-H/{ao) 
XZ 



hrtiao) 



V(a) 



Hrs{a) - Hrs{ao) 1, / X 1, / X 

XZ 2 2 



(4.7) 



where e^fe is the completely antisymmetric unit pseudotensor whose components are zero unless i ^ j ^ k and in this 
case they are equal to 1 or —1 according to the fact that ijk is a even or odd permutation of 123. 

If the electromagnetic field which is under investigation is localized in a space region with linear dimensions small 
compared with 1/x (for instance the electromagnetic field of a resonant cavity), then it is appropriate to expand the 
relations (|4.4| )-(4.7) in power of x^^ ■ The first three terms of the expansion are given by: 



X^x^x^ 
4 



(4.8) 



i^03 = ^03 - ^h^'K' i^Os - ^ss) + 



1 _ 1 



(4.9) 



Fi2 - ^12 - ^h^^'l'e^tsS'' {Tos - ^.3) + 



x'^x^x^ 



6 



-ets3 



(4.10) 
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^^;in.(.^ot-^«)+ (4.11) 



i n wh ich explicit hrs and its derivatives are functions of (— xx") and the tensor part oiTf^, in the first term of (|4.§|) 



( 4.11 ) is to be expanded in power of x until second order. We see immediately that the z components of both 
electric and magnetic fields have the first term of the expansion vanishing. Therefore in the first approximation the 
difference with the TT-gauge arises in the electromagnetic field component normal to the direction of propagation of 
the gravitational wave. 

If the problem is to "follow" a photon in its motion, that is to say to consider the phase shift and the change in the 



polarization of a photon propagating in a gravitational wave field in a region small enough to use Eqs. ( |4.8D -(4.11), 
then for consistency reasons we must also expand those equations in powers of ■ lu this case we have: 

F -T ^h^T 4-^^V/,(i)-F 4-;,(i) T \ (1) 

-f^Ot — -rot — -j^n-t -^Or H Y' \ r -^st + n rt-^OS J ^ " t -^Or 

- ^^h(\iTo, - ^,3) + ^^^^^^h^V^st + (4.12) 
X2a;''(a;3-3a;0),(2) ^ x'(2a;3x3 - Gx^x" + Sx^x") , (2) 

H g 'I rf^03 ri t J-Qr 

Fo3 = ^03 - ^h('->/ {:Fos - :f,ss) ~ ^^h^^l^T,, + 

2 D 



^— -h'^^K'Tos + ^ -h'^^\'Ts3 (4.13) 



r 2 r 3 

- ""'"^^"['^""^ h^Vest^S'^i^o, - T,,) (4.14) 



xVx^,(2) ,^ ^ ^ x^x-(2x^-3x°) (2).^ , 

xVV-3xO) (2) _ x'(2x3x3-6x3x° + 3x°x°) ^(i) 
H ^ rt-^03 n ^ J-r3 



(4.15) 



in which now explicit hrs and its derivatives are intended calculated in zero and also the tensor part of J-'^i, in the 



first term of (4.12)-(4.15) is to be expanded in power of x {x^ ~ 2;°) until second order. 



V. PHASE-SHIFT 

As a first application of the obtained results we calculate the total phase-shift between the light beams in an inter- 
ferometer with arms in the directions defined by the azimuthal and polar angles {9r, (j)r) {r — 1,2). The gravitational 
wave propagates in the positive z direction. The total phase-shift is obtained by the difference of the phase change 
along the two arms: 

= Suif - Su2f + Suib — 6u2b (5.1) 
where lower indices 1, 2 refers to the two arms, and /, b to forward and backward path. 
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The phase is a scalar quantity, therefore we can perform the calculation in TT gauge; the result holds true also in 
FNC, the laboratory frame. In TT gauge, free falling bodies initially at rest, stay at rest at any subsequent time (as it 
is easily seen by inspection of motion equations). Therefore, if the beam splitter is in the origin O of the coordinates 
and the two mirrors are placed at a distance L alon g the directions — (sin 0^ cos 1/)^, sin 0^ sin 0^, cos 0,.) {1' — 1,2), 
the forward phase change is given by (see eq. (2.31)): 



1 k 



^Urf = --(1 + cos6ir) [i?0,((?!)o - - X^(l - cosOr)) ~ H^^{(t>o - xy°)] 



(5.2) 



where i?^ — cos20i7+ + sin20i7x- When light is travelling backward the reverse direction of propagation is given 
by the angles (j: — Oj-^tt + 4'r)- Therefore the backward phase change is given by 



1 k 

Surb = --(1 - cose'r) [H^^{(j)o - X2/° - 2xL) - i^0,(0o - xy" - Xi(l - cos6'r))] 



(5.3) 



This holds true for any waveform of the gravitational radiation. Now we consider a plane monochromatic gravitational 
wave. We set 



^-^(^g) = sin ($g + a); 



= cos 20 j4+ + sin 2(/) 



(5.4) 



where the constants Ax are the polarization amplitudes (/i_|_(<i>g) = sin (<i>g + a), ft-x (,^g) = A^ sin ($g + a)). 
In this case we have 



Surf — ^0^ — (1 + cos6'r) sin 
X 



(1 — cos 9r ) 



sm 



durb = A^i, — (1 — cos6'r) siu sin (xL + xy" — a) 
X 



A^^ — (1 — cos(?r) sin 
X 



xi 



(1 — C0s6'r) 



sm 



xL 



xL 



(1 — cos 9r) + XlP ~ 



(1 — COS 0r) + XiP ~ 



Therefore the phase change in each arm is given by 



Xi 



(1 — COsOr) 



sm 



xL 



6urf + Surb — A^^— <^2 COS 9r sin 

+ (1 — cos^r) sinxi sin (xi + xy*' ^ ct)} 
For small 9r angle, when one arm direction is very near to that of z axis one obtains 



(1 — cos6'r) + xy° ^ 



SUrf + 6Urb 



e,-^o . k 9^ 



- if [sin xL sin {xL + xy° - a) + xL sin (xy" - a)] 
X ^ 



while for almost antiparallel direction one has 



Su. 



rf 



SUr 



A, 



k (tt 



X 



[sin xL sin {xL + xy° — a) + xL sin {2xL + XV^ ~ ") 



We see therefore that for parallel and antiparallel directions there is no phase change. 
In the limit of small distances x-^ << 1 one has 

5urj + 6urb ■^^^^^ A^^ k L SIT? 9 sin (xy° — a) 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



Now let us consider an interferometer with arms in the positive x and y directions respectively. In this case 
Aff,^ = Aj^ and = — A+. Therefore the total phase-shift is given by 



(5$ = 2{5uif + Suib) = A+k (2L) 



sinxi 
xL 



sin {xL + xy° - a) 



(5.9) 



which is in agreement with known results (see for instance p8[ 
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VI. POLARIZATION 



There are two possible state of polarization. If the electromagnetic wave is propagating in the k = 
k{sin9cos(j),sin6sm(f),cos9) direction we are naturally led to describe the electromagnetic field in a frame whose 
three axis directions are: 



6(1) = (— cos 9 cos (f), — cos 9 sin cj), sin 9) 
6(2) = (sin (/), — cos (f), 0) 
fe 

6(3) = T = (sin ^ COS 0, sin ^ sin ^, COS ^) 



(6.1) 



In flat space-time, the possible polarization of the electric field are along the e(i) and 6(2) directions. They are 
obtained when 



■ COS( 



sm( 



El 
k ' 



B^i) = — sin( 

S(2) = COS (/> — 



(6.2) 

(6.3) 



When a polarized wave interacts with a gravitational wave, the electric field rotates and changes in magnitude. We 

can write 



{Tko)^r) = {E^)(rA^) + (^2)(.)e^2) + {E:,\,)e\z) ns + u) 



(6.4) 



in which r = 1,2 and 



1+2 ('^°^ ^ ' 



(-£'i)(i) = 
(£^3)(i) = ^sm9{l-cose)h^ 



' -l)h^- — (cos (j)bl + sin (j)bl) 
-tin 



1 k 

2^0 + (sin#J - cos (f>bl) 



(6.5) 



(£'i)(2) = -Bo 



{E2){2) — Eo 

{E^)i2) = -Eosm9h'^ 



1 k 

- (1 + 2 cos 9)h'^ - — (cos #0 + sin #^) 

1 k 
1 - - (2 + cos 6*) /i^ - — (sin </>' ' 



COS( 



(6.6) 



where 



= cos 2(j) h+ + sin 2^ hx 
h'j, = — sin2(^ /i+ + cos2(^ /ix 



(6.7) 
(6.8) 



As an example of application, let us consider the following "gedanken Experiment". A system is mack; by a 
polarizer in the origin of a TT reference frame and a particle with mass m and charge g at a distance L along 
the direction {9,(f)): both of them are freely falling and point-like with respect to the gravitational wavelength. An 
electromagnetic wave pulse is propagating in the same direction. The duration Tp of the pulse fulfills the condition 
l/{ck) << Tp << l/(cx), namely it is very small compared with the gravitational wave period but it is very long with 
respect to the electromagnetic period. In this way the electromagnetic pulse is short enough that the metric tensor 
could be considered constant in the region in which the electromagnetic field exists but also long enough that this 
could be considered practically a plane monochromatic wave. A charged particle in the presence of an electromagnetic 
field in curved space-time undergoes a force whose expression is 



mc 



(6.9) 
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where ds is the Hne element, u^^ the four-velocity of the particle and F'^^ are the ChristofFel symbols. If the particle 
is at rest in TT gauge, then — (1,0, 0, 0). In this way one obtains 



ds2 



= 



dy° = ds 



mc 

q 



-pi 



(6.10) 



m c 



pr 



Let us suppose that the polarization is in the 6(2) direction and that the electromagnetic beam reach the polarizer at 
the time U. Therefore at the particle position 



= --Eo sine [h'^iaf] - h'^{a,)] /'[k^y'^ + u{af)] 



(6.11) 



where a/ = x{L cos 6* — ti — L), Ui = x{Lcos9 — ti) and e [ti + L.ti + L + Tp], e L e^g^, {L + Tp) e^g^ 
particle has therefore a TT acceleration in the direction z given by 



df^ 



2m 



■ sm 



e[h'^{af)-h'4ai)] fik^y'^ + uiaf)] 



The 



(6.12) 



This holds true for any waveform of the gravitational radiation. Now we consider a plane monochromatic gravitational 
wave [sec eqs. (5^)]. Therefore [see also Eq. (6^)] h'A^g) = A'.sii\{^g + a), A'. — — sin2(/) ^+ +cos20 A^, and 



df^ 



QEoA' 



■ sm W sm 



^(1-COS0) 



cos 



a - xti - ^(1 - cost 



(6.13) 



It is immediately seen that for parallel and antiparallel interaction the acceleration vanishes, 
interaction one has 



'2„,3 



For xL « 1 



dy_ 

df^ 



df^ 



qEoA' 



qEpA'^xL 
2m 



xL 
2 



cos 



Xt. + ^-a^/'[fcpy^ + "(a/)] 



For perpendicular 



(6.14) 



sin0(l — COS0) cos [xti — q;)/'[/c^j/^ + u(a/)] 



(6.15) 



from which one could easily see that the greater acceleration is caused when = 2 7r/3. In order to get the observable 
acceleration, one must get its expression in FNC, being this the laboratory reference frame. However the acceleration 
is a first order quantity in /i^^; therefore its expression in FNC is the same as in TT gauge. 

It has been shown that, in the laboratory reference frame, the particle is accelerated in the z direction only because 
of the direct interaction between the incoming electromagnetic and gravitational waves. In the absence of one of 
the two wave the acceleration would vanish, thus being a peculiar feature of the electromagnetic propagation in a 
gravitational background. This result contains an apparent paradox: the appearance of a component of the electric 
field along the gravitational wave propagation direction with the consequent acceleration of the particle along the 
wave vector of the gravitational wave, even though this last one is a transverse wave. This is due to the fact that the 
components of the electromagnetic tensors are not independent (because of M a xwell equations). 



In order to get the order of magnitude of the effect described by Eqs. ( |6.12 )-( |6!T5| ) let us suppose that the charged 
particle is a proton. Taking = 7r/2, a laser beam with intensity / = 50 W/mm^ and wavelength A = 271 /k ~ 1 /im, 
and putting the particle at a distance L = t:/x from the polarizer (in such a way to have the greatest acceleration) 
we find that, after the passage of the electromagnetic wave pulse, the drift velocity in the z direction is [the factor 7 
takes into account the effect of the phases of the last two factors in Eq. ( 6.13| )] 



= 7 0.7 A'm cm/sec. 



-i<7<i; 



(6.16) 



thus, even in this optimistic situation, the effect is exceedingly small and, most likely, out of the range of an experi- 
mental verification. 
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As a second example of application let us consider the possibility to perform a polarization rotation measure. 
PVLAS is an experiment under construction which is designed to measure the vacuum magnetic birefringence 



The achievable sensitivity of about 5 10^^ rad/\/Hz, with an integration time of 10^ sec |30 , would allow the 
measurement of an ellipticity of 10~^^ rad. With a simple modification of the experimental setup (the collocation of a 
quarter-wave plate) the same sensitivity could be reached in the rotation angle of a linearly polarized electromagnetic 



radiation |30|. This is the best result to date. The description of the apparatus is found in |29|; however, in the last 
years there have been many improvements in the performances. For instance the finesse of the 4.5 m Fabry-Perot 
cavity is expected to be 10^ 10^, while the laser have an intensity of ~ 100 mW/mm^ In the same assumptions 
as in the previous example for the ele ctrom agnetic pulse, and for gravitational waves of ~ 100 Hz, the approximation 



xL <C 1 could be used, and then Eq. ( 6.11 ) gives on the analyzer prism 



Fo3 ^Toz^F^E^ v4; sine cos {xU - a)f[k^y^' + u{af)] (6.17) 

where F is the finesse of the cavity and the first equality holds for first order quantities in hf^i, . The rotation angle is 
therefore given by 

sine F ^ A'^^ sine A'^. (6.18) 

It is therefore seen once again that also taking the most optimistic value of 10^^"* for A'^ (gravitational wave from a 
pulsar) the effect is too small to be detected. 

By the way we note that in Ref. ||2^] authors claimed that no rotation of the plane of the polarization of the 
electromagnetic wave occurs to first order in h calculations. Indeed this result was obtained for a particular linear 
polarization of the gravitational wave. In fact authors set = 0. In this particular case also our solution is 
consistent with the previous claim. However in the general case there is a rotation of the plane of the polarization of 
the electromagnetic wave. A comparison showing more in detail the differences between the approach of this paper 
and the splitting procedure in the physical situation described in ||2^ is found in 



VII. CONCLUSION 



We have found a solution in TT gauge to the free de Rahm equations with Lorentz condition for a gravitational 
background described by a plane wave, within the framework of two approximations. The first one is the linear gravi- 
tational approximation. The second one is that the gravitational wavelength is much smaller than the electromagnetic 
one; this allowed us to find a solution one order of approximation beyond the limit of geometrical optics. The result 
was used to write the components of the electromagnetic tensor in TT gauge and FNC which automatically fulfill free 
Maxwell equations in curved space-time up to terms of order {x/k)^ . 

We have applied these solutions to the problem of polarization of an electromagnetic field, exploiting an interesting 
feature of the propagation of electromagnetic waves in a gravitational wave background. In fact we have found that, 
when an electromagnetic field linearly polarized interacts with a gravitational wave whose direction of propagation is 
perpendicular to the electric field, then the electromagnetic field changes in magnitude and rotates in such a way that 
the electric field gets a component parallel to the gravitational wave vector. This result could seem in disagreement 
with the transversal nature of the gravitational waves: however this fact arise from the mutual dependence of the 
electromagnetic tensor components through Maxwell equations. 

Two possible applications of this result have been investigated. However the effects are so small that are out of the 
range of current experimental techniques, and unlikely to be ever observed in the future. We think that a study is 
recommended in order to check if possible measurable effects could be found on light beam coming from cosmological 
distances (analogously to what was done in fl^ for amplitude fluctuations). 
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